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Abstract - In this work, we present the analytical solution of the effective mass Pauli equation, 
with Rashba and linear Dresselhaus interactions, for an electron gas moving through a semicon- 
ductor quantum dot under a longitudinal electric field, which is defined along the a;-direction. We 
study the relative influence of the Rashba and Dresselhaus terms on the spin-Hall effect for the 
first propagating and edge channels, by analyzing the mixing between spin-up and -down states 
and the zero-field spin splitting along the transverse directions. When the spin rotation depends 
only on the {/-coordinate, the spin orientation and the spin density vary along this transverse 
coordinate and, in this case, we show that the spin-Hall effect is only due to the Dresselhaus term, 
for depolarized electrons. On the other hand, if the spin rotation depends on the z-coordinate, 
the spin-Hall effect is provoked only by the Rashba interaction. 



Introduction. — In a semiconductor, a two- 
dimensional electron gas (2DEG) in the conduction band 
does not experience as strong a nuclear attraction as an 
electron in an atom does. However, it may still see electric 
fields due to internal effects. Those symmetry-breaking 
electric fields, on one hand, can induce a spin splitting 
due to a structural inversion asymmetry, which is named 
Rashba spin-orbit interaction [UIl]- On the other hand, 
the well-known Dresselhaus spin-orbit interaction [3] can 
induce a splitting by electric fields that arise due to bulk 
(or crystallographic) inversion asymmetry. Externally ap- 
plied electric fields may also provoke such spin splittings 
[^{8]. The full understanding of the spin dynamics in 
a semiconductor heterostructure is of great interest from 
the theoretical and the applied physics points of view [9] . 
Taking into account both spin-orbit interactions (SOI), 
theoretical work has been focused on the investigation of 
electron transport in a 2DEG j¥^ET|. We propose here 
to solve analytically the quasi two-dimensional Pauli eq. 
with a longitudinal electric field, which induces both SOI. 

Among other experiments, in a 2DEG in strong mag- 
netic fields, the integer Quantum Hall effect and the Hall 
resistance quantization [52] lead to transverse modes and 
edge states. In this case, each edge state is spin-degenerate 
because the Zeeman splitting is negligible [23l[24]. Re- 



cently, theoreticians predict a different kind of edge chan- 
nels that can form even in the absence of an external mag- 
netic field j4l[5l [25H30] ■ this anomalous behavior is called 
spin-Hall effect (SHE). There are experimental demonstra- 
tions of the SHE in quantum wells of HgTe [31] and of 
InAs/GaSb [32]. 

In this work, we study the dynamics of electrons in 
quasi two-dimensional quantum dots in the presence of 
Rashba and Dresselhaus spin-orbit interactions, which are 
provoked by an applied electric field. In this approach 
(see the next section), we neglect the cubic Dresselhaus 
spin-orbit term. Here, we work with regions where there 
is no applied magnetic field. We place special emphasis 
on how to use the "variable separation method" since we 
have written the eigenspinors as the product of a trans- 
verse space dependent part (a particle eigenfunction in an 
infinite potential well) and a spin dependent part (a ro- 
tation operator times another spinor) for each transverse 
direction. The solution of the longitudinal part is given 
by the well-known Airy functions. We find spin accumula- 
tions near the semiconductor quantum dot edges: spins of 
one polarization pile up at one edge, while the spins with 
opposite polarization pile up at the other edge. In other 
words, it is possible to predict the spin edge states and, 
therefore, to give a description of SHE. 
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The Pauli equation and its solutions. To study 
transport properties of a quasi-2DEG in a longitudinal 
electric field, we shall consider a quasi two-dimensional 
quantum dot of dimensions > Wy > w^, such that 

2 2 2 

< ^ < We will solve the effective mass Pauli 
eq. for an electron moving in the direction x under a 
longitudinal electric field F = 



+Hr + Hd 



eTx + V{x, y, z) 
- Ey] * {x,y,z) : 



0, 



(1) 



where # {x, y, z) is a spinor, p = —ih^ is the momen- 
tum, m is the effective mass, is the Fermi energy, 
V{x, y, z) is the confining hard wall potential and Hr and 
Hd are the Rashba and the Dresselhaus interactions re- 
spectively. We consider that both SOI arc induced by the 
applied electric field J- and, therefore, Hr has the follow- 
ing form Hr = J- {(JyPz — (^zPy), with a* ~ aeT and a 
the material constant. On the other hand, it is possible 
to show that the average momentum squared in ground 
state is given by (p^) = 0.7794 {Ti/lEf, where 1% = 
is the electric length Notice that for high values of 
the electric field, such that > Ie, the following in- 
equalities (p^) > (p^) > (Py) can be true and we could 
propose that the most important term of the linear Dres- 
selhaus interaction is written by Hd = {cTyPy — ctzPz), 
with P* ~ 0.7794:f33D/lE a-nd /Jsd is the Dresselhaus three 
dimensional term. 

The Pauli eq. under our approach is reduced to the 
following form 



+2i [aay - Pal) ^] * (.t, y, z) = 0, 



(2) 



here a = ma* /h?' and /3 = raf3* /h?' . In order to get rid of 
the variable 2, we will consider the complete set of station- 
ary states of a particle in a z-dimensional infinite potential 
well with boundary conditions Z„(— = Z^iw^l'^) = 
0, which are solutions of the differential eq. 



It is easy to verify that C En d = En (^) ^n, 
where C,n and c„ are spinors. 



(z) 



cos (fc^.n z) for n odd 
sin (fc^^n z) for n even 



h — nJL n 



-a'^+P'^ and = e 



-i{ai7y — l3(T~)z 



rotation operator. If we use the complete set of functions 
{Cn} to expand $ (x, y, z), we get 



* (x, y, z) = Er=i Zr,{z)Fn{x, y). 



(3) 



By introducing this function in the Pauli eq. ([2|), multi- 
plying by R|Zi (z) and integrating on the variable z, we 



have 



2mE-p 



7f 



-^(2/32 + (a2-/?2)JM)]^}F, 
;y2^ En [W^^n,/ - W'^y - acTz) ^,/] ^ = 



(4) 



here rj = ^ o? Ar and the mixing terms are given 



by In.l 
2 



■r 



-./2 



' ^2 - i.2 



l^sa\2r]z Zn{z)Zi(z)dz and J„_; = 

J^^^^2 cos 277Z Zn{z)Zi{z)dz. Notice that if /3 = ±a 
the previous eqs. are uncoupled. 

When the width Wz is narrow, the energy levels spac- 
ing, that is defined by lS.Ez,i — 

{21 + 1), is wide. Therefore, the transverse reso- 
nances are clearly separated in energy in the quasi two- 
dimensional quantum dot and rjz < rjWz << 1. We 
are able to take the approximations sin 2772 ~ 2riz and 
COS277Z ~ 1, in such case /14 ~ 0, Jia — 1, h,2 — ^g^^' 
and J1.2 — 0. In the case where there is only one propa- 
gating mode (for Fermi energies Ep < Ai?i), we are left 
with the 2D differential eq. 



'-Ef 



+2i {P<jy -aaz)^ 



wt 
dy 



= 0, 



(5) 



which describes the quantum dot in the 2DEG limit. We 

« 1. 



"9^ 



neglect here the term /i^i, because 
For clarity, the index n = 1 shall be suppressed in the fol- 
lowing expressions. Taking into account the ansatz given 
by ([3]), we use a similar one 



F{x,y) = RySYZiYiiv)M^ 



(6) 



where R,, 



^-i(l3<yy-a!7^)y jg auothcr rotation operator 
acting on the yz plane, Xi is a scalar function, Yi{y) is 
an cigenstate of a particle in an infinite potential well 
and S' is a spinor, which is independent of x. The 
sum runs over the N propagating modes, with trans- 
verse wave number ; = (iTr/wy) + (ir/wz) such that 
n'ff = — K± I > 0, and an infinite number of evanes- 

cent modes. For each longitudinal wavcnumber ni there 
are two propagating physical channels: one with spin- up 
and another with spin-down. Equation ([5]) gets trans- 
formed into the following system of uncoupled eqs. 



2if 



(7) 



for each Fourier open channel / at the y direction. The 
general solutions of these eqs. are given by 



Xi (x) = Ai(-Xi)Co+Bi(-Xi)Ci, 



being 



Ai(-xO 



/X; 
3 



(|xf)+Ji(|x? 
-/-Jfx?)-^! fx? 



(8) 



(9) 



p-2 
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the well-know Airy functions and = + + 2?^^) l^- 
If we follow the procedure developed in [5], the solution 
([8]) is reduced and it has the following form 



Xi^i (x) = 1A2611 



-1/2 



ri 



with 7'i = —2.3381 as the first zero of Ai and E^.i = 
—rieTlE as the ground state energy. Here, the Fermi en- 
ergy is totally quantized fi'F.i,; = + i- From now 
on, we will suppress the index 1 for all expressions. Notice 
that those eqs. indicate that there is no mixing among the 
^-Fourier channels. Due to the rotation operator Ry, there 
is mixing between the spin channels. 

Edge states as the origin of the spin-Hall effect. 

— The spin-Hall effect is due to spin accumulations near 
the semiconductor quantum-dot edges even in the absence 
of an external magnetic field. In this section, we will show 
that depending on the Dresselhaus interaction, these rota- 
tions might provoke spin accumulations close to y = 
or to y = for electronic depolarized currents. We shall 
analyze the principal characteristics of the spin-transverse 
wave functions at a particular case where there are depo- 
larized electrons. 

According to the previous section, the spin-dependent 
transverse wave function is written by 



Fiix,y)^RySYiiy)Xi (x) 



(10) 



The exponential term Ry = e~'^^P''y~°"^'>y is a rotation op- 
erator and indicates a precession about an axis in the plane 
yz. Due to the Rashba interaction, electron spins rotate 
around the z axis, while the Dresselhaus term causes rota- 
tion around the y axis. The mixing between spin channels 
depends only on the strength of the Dresselhaus term. 

It is well known that the GaSb has the following char- 
acteristics: a = 33A2, B^D = 187 eVA^ and ^ = 0.0412 
[6]. On the other hand, we use the following param- 
eters: = 5.0 X 10^^, = Wy = 60 nm and 
Wz = 20 nm for the first longitudinal ground-state energy 
Exs = -rieFlE = 0.3345eV. In figs. [T]and[l]we plot 
the spinors' amplitudes |-F/,^|^ and the difference 

between them A \Fi^ = ~ as functions of the 

y direction, for Z = 1 and 2. 

Figure [1] shows the spin amplitudes for spin depolarized 
electrons {Si^-\ = l/V^ and Si^i = l/V^) considering that 
there is no Dresselhaus term. For the sake of clarity, the 
curves for |Fi,tr and \Fi. ^1 are vertically offset 0.1 am- 
plitude units. The electron spins precess only around the 
z direction, preserving their own initial polarization, be- 
cause there is only the Rashba interaction, thus one can 
not find any spin accumulation near the edges. In this 
way, the Rashba term does not affect the symmetry of 



iFi^^i" and 



Yi^ (y) for 



because \Fi^^f = |^/,;|^ 
I — 1,2. Therefore, there is no difference between the 
spinor amplitudes: on both cases AlFil =0. 

Let us now cancel the Rashba contribution. Figure [2] 
shows that for spin depolarized electrons the symmetry of 
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Fig. 1: l-Fi.tl , andA|_Fi| as functions of the j/ direction 

when there is no Dresselhaus term /3 = and a 7^ 0, for spin 
depolarized electrons at the 2 direction: a) when I = \\ and 
b) when I = 2. For the sake of clarity, the curves for 
and l-Fi.^l^ are vertically offset by 0.1 amplitude units. The 
electron spins only precede about z direction, preserving their 
own initial polarization, and this motivates that one can not 
find any spin accumulation near the edges. 



and arc broken by the Dresselhaus term. In 

this case, the exponential term is reduced to the following 
form 



(11) 



g-m„,9y ^ COS {I3y) -sin(/?y) 
|_ sin (/??/) cos(/3j/) 

this exponential term indicates mixing between both spin 
polarizations because the off-diagonal terms are real and 
different from zero for processes where a spin flip has taken 
place in the transverse y-direction; and, therefore, the spin 
amplitudes are given by 



i [1- Sin (2/32/)] 



[1 



(y) 

sin (2/32/)] (y) 



where 



Yi{y) 



cos (kyjy) 
sm{ky^iy) 



for I odd 
for / even 



(12) 



(13) 



and kyj = The factor corresponding to a stationary 

particle in a infinite potential well YJ^ (y) dominates at 
the edges' neighborhood and at its possible nodes, caus- 
ing suppression of the spin-dependent amplitudes in such 

I 1 2 

zones. In the neighborhood of 2/ = 0, are dominated 



p-3 
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Fig. 2: \Fi^ff , |-F1,4.|^ and A|f;|^ as functions of the y direc- 
tion when here is no Rashba term, a = and /3 7^ for spin 
depolarized electrons at the z-direction: a) when I — 1; and b) 
when I = 2. The symmetry of and 1^1,4.]^ is broken by 

the Dresselhaus term. The spin-up electrons pile up at the left 
edge because A > 0, while A < indicating that the 
spin-down electrons pile up at the right edge. 



by the 1 ± sin {2/3y) factors. Both factors are positive and 
1 — sin (2/3j/) increases while 1 -t- sin {2/3y) decreases when 
y grows. The maxima of |i^(,-t-| are moved and grow to 
the left side by 1 — sin {2/3y), while for maxima are 

moved to the right side due to the 1 + sin (2/3y) term. Wc 
must notice, however, that A|F;|^ = — sin (2/3y) YJ^ (y) 
and, therefore, the spin-up electrons pile up at the left 
edge because A \Fi\^ > 0, while A \Fi\^ < indicates that 
the spin-down electrons pile up at the right edge, these 
spin accumulations imply a mixing between the spin pop- 
ulations. In such cases, \Fi^^f and l-F;,^]^ become the well- 
known edge states. 

The behavior is described in figs. [T] and [5J shown only 
for the / = 1, 2 states. For other states the general trends 
are conserved. In this way, wc only discussed them for 
; = 1 and 2. 

When there are both SOI, the exponential term now is 
written in the following form 



-i{Pay-aij^)y _ 

COS (?7y) sin {T]y) 
T: sin ivy) 



sin (to) 
cos (to) - i^sin (to) 



(14) 



the off-diagonal matrix elements sin {rjy) reflect the 
passage of flux from one spin state to another. The ori- 
gin of these transitions is the Dresselhaus spin-orbit inter- 
action that stimulates a mixing between the propagating 
physical channels, in the transverse regions where the /? is 
significant. On the other hand, the spin-dependent edge 



states are given by = 
^ sin (2to) 



1 



sin (2?/?/) Yj^ (y) and 



F, 



1 



iy). 



The eigen- energies, the spin-polarized wave functions 
and the zero-field spin splitting on the y-direction can be 
obtained even if we restrict our analysis to energies bellow 
the N propagating mode threshold Ef < ti^ Kjf_^_^/2m = 

^ {{N + l)TT/wy)^ + (tt/wz)'^ and neglect the evanes- 
cent mode contributions. In this way, we found the an- 
alytic solutions given by eq. pH)) . This solution refiects 
the mixing of spin f and spin J, states in the j/z-plane. 
The propagating modes of the spin polarized states are 
obtained after a unitary transformation 



U 



Here, 



Fi-^ Fi = Ue-''^'^''y-"'^-'>yU''USYiXi (x) 
Fi{x,y)^e"'''>ySYi{y)Xi{x). 



(15) 



u = ut 



[{a + r])az- fidy] 



is the relationship of the similarity matrix, Fi(x,y) = 
\JFi{x,y) and S = US. In this representation for the 
right-side propagating state spinors at any point y in- 
side the wave guide, the energies of the spin f and spin 
I states are E-^ ^ [^^.r + — 2ky,j.'>l\ ^i^d E\^ = 
[^y r + '7^ + 2ky^r'']] ^^'^ thc zero-field spin splitting is 
given by E^ - E^ = ^ky^rV- 

On the other hand, taking into account the 2D effective 
mass Pauli cq. ([5]) we can define an electric current density 
j that guarantees the conservation of probability. The y 
component of such current is given by 



■ _ -leh 
■ly^l- ~ 2m 



dy 



!pf dFi 
■I dy 



2iF{ {p(jy - an,) Fi 



with 



Fi{x,y) - RySYi{y)Xi{x) 



(16) 



the complete analytic solution. Notice that jy^i = for 
Fi. This means that, at variance with the classical Hall 
effect, no Hall voltage builds up between the two edges 
since no charge imbalance results from this phenomena at 
the y-transverse direction. This spin imbalance is part of 
the spin texturing and forms an anomalous Hall effect, the 
spin-Hall effect. However, given that quantum dot is made 
of a non-magnetic semiconductor, it has equal populations 
of electrons with spin-f and -J,. 

To solve eq. ^ we could follow a different procedure 
that leads to a differential eq. in only one variable x similar 
to d?]). Let us suppose that the quasi-2DEG is in the xz- 
plane. Now Wx ^ Wz > Wy. To get rid of the variable y. 



p-4 
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we will consider the complete set of stationary states of 
a particle in a y-dimensional infinite potential well and, 
following similar steps to get ([5]), we obtain 



Gi 



0, 



(17) 



Taking into account the ansatz given by we use a 
similar one 



G{x,z) 



SYZiZi{z)Xi{x), (18) 



and so on. It is clear that, in this case, the Rashba inter- 
action provokes edge states at the z-direction. 

Summary and conclusions. We study the trans- 
port properties of spin-^ electrons in an homogeneous 
semiconductor quantum dot subject to an externally longi- 
tudinal electric field. By describing the interacting evolu- 
tion of these electrons with Rashba and linear Dresselhaus 
interactions, we found the analytic solution of the Pauli 
eq. for spin-up and spin-down coupled channels and de- 
termined the spin texturing; in other words, the spin ori- 
entation and the spin density can vary along the transver- 
sal directions. On the basis of this phenomena, we found 
the edge states. Such solutions arc in terms of a rotation 
operator and it indicated a precession about an axis in 
the plane j/z, rotations that might provoke spin accumu- 
lations close to the edges. There are no charge imbalance 
resulting from this phenomena at the y-transverse direc- 
tion, because the transverse current is zero and, therefore, 
there is no Hall voltage build-up between the two edges. 

We can propose, alternatively, that the slab which de- 
fines the quantum dot is in the xz-plane, here > >> 
Wy. The 2D effective mass Pauli eq. now is in terms of the 
transverse z-direction. When we solve it, taking into ac- 
count the rotation operator e"*'^"'^" ~^'^»=)^, the edge states 
and the mixing between both spin polarizations are pro- 
voked only by the Rashba term. 
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